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Abstract 

THIS PAPER IS ELIGIBLE FOR THE STUDENT BEST PAPER AWARD. We consider 
in-network computation of MAX in a structure-free random multihop wireless network. Nodes 
do not know their relative or absolute locations and use the Aloha MAC protocol. For one-shot 
computation, we describe a protocol in which the MAX value becomes available at the origin 
in 0{^Jn/ \ogn) slots with high probability. This is within a constant factor of that required 
by the best coordinated protocol. A minimal structure (knowledge of hop-distance from the 
sink) is imposed on the network and with this structure, we describe a protocol for pipelined 
computation of MAX that achieves a rate of ri(l/(log^ n)). 

1 Introduction 

Early work on computation of functions of binary data over wireless networks focused on computing 
over noisy, time-slotted, broadcast networks, e.g., (HE]. With increasing interest in wireless sensor 
networks, recent research has concentrated on 'in-network' computation over multihop wireless 
networks, e.g., O [Sj H], which also assume that time is slotted. The primary focus of the above 
research has been to define an oblivious protocol that identifies the nodes that are to transmit in 
each slot. This implies that the nodes have organized themselves into a network and have their 
clocks synchronized. Both of these require significant effort. In this paper we describe a protocol 
for in-network computation of MAX in a structure-free network, (i.e., in a network where nodes do 
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not have an identity and hence do not know the topology) that uses the Aloha MAC protocol. We 
first describe the One-Shot MAX protocol for one-shot computation of the MAX and its analysis. We 
show that, with high probability (w.h.p.), the sink will have the result in a time that is within a 
constant factor of that required by a structured network. We then impose a minimal structure and 
describe the Pipelined MAX protocol and its analysis. We show that the rate of computing the MAX 
in this network is \ ). 

^ log^ n ' 



2 MAX in Multihop Aloha 

n nodes are uniformly distributed in [0, 1]^ and each node is assumed to know n. The sink, the 
node that is to have the value of the MAX, is at the origin. The nodes do not have an identity and 
they do not know either their relative or their absolute positions. Hence, the network does not 
know its topology. This of course means that a schedule for transmissions cannot be defined. Thus 
a random access protocol is an obvious choice at the MAC layer. We first assume that the nodes 
use the s- Aloha MAC protocol. For pedagogical convenience, we will assume slotted- Aloha at the 
MAC layer. The analysis easily extends to the case of pure Aloha MAC. 

Spatial reuse is analyzed using the well-known protocol model of interference [3]. For s- Aloha, this 
model translates to the following. Consider a transmitter at location x\ transmitting in a slot t. 
A receiver at location X2, can successfully decode this transmission if and only if the following two 
conditions are satisfied. (1) ||x2 — xi|| < r„, and (2) ||x2 — xsH > (1 -|- A')r„ for some constant 
A' > 0; 3:3 is the location of any other node transmitting in slot t. r„ is called the transmission 
radius. A transmission in slot t is deemed successful if all nodes within r„ of the transmitter receive 
it without collision. The following is a sufficient condition for successful transmission by a node 
located at x in a slot: ||x — > (1 -|- A)r„, A = 1 -|- A', for every other node transmitting in that 
slot and located at x' . 



2.1 One-shot computation of MAX using Aloha 

Let Zi be the value of the one-bit data at Node i and Z := maxi<j<„ Zj. The protocol One-Shot 
MAX is as follows. Node i can either receive or transmit in a slot but not both. In slot t, Node i 
will either transmit, with probability p or listen, with probability (1 — p), independently of all 
the other transmissions in the network. Let be the value of the bit received (i.e., correctly 

decoded in the absence of a collision) by Node i in slot t, t = 1,2, ... . If Node i transmits in 
slot t or if it senses a collision or idle in the slot, then it sets Xj(t) = 0. Define 1^(0) = Zi and 
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Yi{t) := max{li(t - l),Xi{t)} for t = 1,2, . . . . Yi{t) is the 'running MAX' at Node i in slot t. If 
Node i transmits in slot t, it will transmit Tj(t) = — 1). 

It is easy to see that the correct value of Z will 'diffuse' in the network in every slot. The performance 
of the protocol, that is, the diffusion time, depends on p. The choice of p is discussed in Section |3| 

To study the progress of the diffusion, we will consider a tessellation of the unit square into square 
cells of side s„ = r^ 2.75'iogn l '^^^^ ^^^^ result in /„ := ^ = [^/ ^ y^"^^^ -] rows (and columns) 
of cells in [0, 1]^. There will be a total of M„ := \ = / „ 7c?„p.„ l ^ cells. Let C denote the set of 
cells under this tessellation. Let Sc be the set of nodes in Cell c and Nc be the number of nodes in 
Cell c. Under this tessellation, two cells are said to be adjacent if they have a common edge. Let 
the transmission radius be r„ = ^Jl^^l^^El}: ~ ^/bsn- For this value of r„ the network is connected 
w.h.p. [3 . The expected number of nodes in a cell is ns^ ~ 2.75 log n. Further, from Lemma 3.1 of 
[T], for our choice of r„ and Sn, 

Pr (ci log n < A^c < C2 log n for 1 < c < M„) 1 (1) 

where ci = 0.091 and C2 = 5.41. Our results will hold for networks that are connected and which 
satisfy ([T]). From the choice of r„ (i.e., r„ > \/5sn) a successful transmission by any node from 
Cell c is correctly decoded by all nodes in Cell c as well as by all nodes in cells adjacent to Cell c. 

The value of Z can reach the sink along any of the many possible trees rooted at the sink. For our 
analysis, we will divide the progress of the diffusion into the following three 'phases' and analyze 
each of the three phases separately. We reiterate that the above sequence of phases is for the 
purpose of analysis of the time to diffusion. The nodes do not perform any such organization. 



• Phase I for data aggregation within each cell. This phase is completed when every node of 
the network has transmitted successfully at least once. 

• Phase II for progress to the bottom of the square. In this phase, the locally computed values 
of the MAX get diffused into the cells on one side of the unit square as shown in Fig. [T] 

• Phase III for progress into the sink. In this phase, the value of MAX reaches the sink at the 
origin in the manner shown in Fig. [T] 

We show in Section [s] that Phase I will be completed in 0(log^ n) slots w.h.p.. Phase II and Phase 
III will each be completed in O ^y^'logn) slots w.h.p. These results are combined into the following 
theorem. 
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Sink located at the origin 



^ Phase II *■ Phase III 

Figure 1: Direction of diffusion during Phase II and Phase III of protocol One-Shot MAX 

Theorem 1 // all the nodes execute the protocol One-Shot MAX, then for any a, A; > 0, the maxi- 
mum of the binary data at the n nodes is available at the sink with probability at least (l — ^) in 
O (v^) ^lots. 



We will also argue that in another O number of slots, the value of Z would have diffused 

to each node of the network. Note that the best one-shot protocol in an organized network, under 
this choice of r„ will also require ^^^/j^^^ time slots for a one-shot computation of MAX. The 
bound on the time in Theorem [T] is therefore tight. 



2.2 Pipelined computation of MAX using Aloha 

If Z were to be computed continuously using the One-Shot MAX protocol, a throughput of 

can be achieved. We believe some structure in the network is necessary to do better. We will as- 
sume that all nodes have a transmission range that is exactly r„. This strict requirement can be 
easily relaxed but we will keep this assumption for pedagogical convenience. 

We impose the following structure in the network. Prior to the computation, each node obtains its 
minimum hop distance to the sink. Henceforth, we will refer to this as simply the hop distance of 
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the node. From ([T]), each cell in the tessellation is occupied. Since nodes in adjacent cells differ in 
their hop distance by atmost 1, the largest hop distance of a node in the network is no more than 



Let hi be the hop distance of Node i. Observe that a transmission by Node i can be decoded 
successfully by Node j only 'd \hi — hj\ < 1. Hence, if there is a reception by Node i in slot t, 
then that transmission must have been made by a node with hop distance either {hi — 1), hi, or 
{hi + 1). Thus, if a node transmits its hop distance modulo 3 along with its transmitted bit, then 
every receiver that can decode this transmission successfully, can also, by the receiver's knowledge 
of its own hop distance, correctly identify the hop distance of the transmitter. 

Time is divided into rounds, where each round consists of r slots. Minimizing r would maximize 
the throughput. We will discuss this in Section [3] Data arrives at each node at the beginning of 
each round, that is, at the rate of 1 data bit per round. Let the value of the bit at Node i in the 
round r be Zi{r). Z{r) := maxi<j<„ Zj(r), for r = 1, 2, . . . , is to be made available at the sink 
node. Node s. 

Pipelined MAX protocol is as follows. The sink only receives data and does not transmit. The other 
nodes in the network perform the following. (We remind here that the naming of the nodes is for 
our convenience. The nodes themselves do not know their identity.) 

In each slot. Node i either transmits with probability p or listens with probability {1 — p) indepen- 
dently of all other transmissions in the network. The value of p is chosen as in the One-Shot MAX 
protocol. Each node executes the following protocol for round r. 

Transmission: If Node i transmits in slot t of round r, then it transmits three bits {Ai, Bi,Ti{r, t)) 
in the slot. Bits Ai and Bi are the identification hits and are obtained as {hi mod 3). The bit Tj(r, t) 
is the transmitted data hit and is obtained as 



Here, by convention, Zi{v) = Yi{v) = for v <Q. Yi{r — 1) is computed from succesful receptions 
in round (r — 1), as described below. 

Reception: In round r. Node i maintains Yi{r, t) for t = 0, 1, 2, . . . , r. Yi{r, 0) is initialized to at 
the beginning of round r. Yi{r,t) stores the MAX of the data bits that Node i has decoded from 
all the slots in round r, upto and including slot t, and which were transmitted by the nodes with 
hop distance {hi + 1). In slot t of round r, if Node i successfully receives a transmission from a 
node with hop distance {hi + 1) (available from the identification bits), then it sets the received 
data bit to be Xi{r, t). If Node i senses an idle or a collision in slot t, or if it receives a successful 




Ti{r,t) 



max{Zj(r -d + hi),Yi{r - 1)}. 
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transmission from a node with hop distance different from {hi + 1), then it sets Xi{r,t) = 0. Thus, 
Yi{r, t) = max{yi(r, t - 1), Xi(r, t)}. Define Yi{r) := Yi{r, r). 

The sink node, Node s, obtains the MAX as Z{r — d) = max{Zs(r — d),Ys{r)}^ for all r > d. The 
delay of the protocol is d rounds or dr time slots. 

Theorem 2 // all the nodes execute the protocol Pipelined MAX, then for any a,k > 0, there exists 
T = T{a,k) = 0(log^ n) so that the correct MAX is available at the sink in a round with probability 
atleast (l — ^) . This achieves a throughput of ^ ^ ^ with a delay of n log'^ n) slots. 

The optimal pipelined protocol for MAX in an organized network requires 0(logn) slots for each 
round in the absence of block coding. Thus, the penalty for minimal organization and no coordina- 
tion is the logn overhead for the length of each round. Also, for our protocol. Node with a hop 
distance of hi, requires a memory of {d — hi + 1) bits to store Zi{r), Zi{r — 1), ... , Zi{r — d + hi). 
Thus, the protocol requires each node to have (d + l) bits of memory for storage of past data values. 



3.1 Preliminaries 

3.1.1 Bounding the Number of Interfering Neighbors 

Define the interfering neig hborhood of Node i by := {i : < \\Xi - Xj\\ < (1 + A)r„}. As 
discussed earlier, a transmission from Node i in slot t is deemed successful if all nodes within r„ 
of Node i can decode this transmission without a collision. A sufficient condition for Node i to be 
successful in transmitting in slot t is that no node belonging to M- ^^ must transmit in slot t. 

From the protocol model, the choice of s„ and ([T|, the set of nodes that interfere with a transmission 
from a node in Cell c, (i.e., UiG^c -^i^^) contained within an interference square centered at Cell c. 



3 Proofs 



This square contains ki 




(1+A)r, 



n 



] + 1) cells. From Q 



\K \ ^ kiC2logn-l 



(2) 



Observe that ki is a constant for large enough n. 



Accepted at ISIT 2008 



6 



3.1.2 Probability of a successful transmission from a cell 



Let Pi be the probability that Node i transmits successfully in a slot and P^'^\ the probability that 
some node in Cell c transmits successfully in a slot. Pi > p{l I, and from we have 

Pi ^ p(l — j))'^i^2 iogra-i_ Successful trausmissious by nodes from Cell c are mutually disjoint events, 
and hence, p(^) = Eie5c ^« - ^^^(1 - p)fciC2 iogn.-i_ pxom we have Nc > cilogn Vc G C and 
hence, P*^*^-* > cilogn p{l — p)'=iC2iogr!,-i^ Choosing p = j^^^^i^^^ maximises the lower bound in this 
inequality and yields 

/ 1 \ -(fclC2 logn-l) 

fciC2 V fciC2logn— 1/ A:iC2e 



Thus, the probability of successful transmission from a cell is lower bounded by a constant ps, 
independent of the number of nodes in the network. This will be crucial to our analysis. 



3.2 Proof of Theorem [T] 

We will prove Theorem [T] by proving bounds on the total time required by each of phases I, II and 
III. 



3.2.1 Phase I: Data aggregation within each cell 

Consider Cell c. Let Tc be the total number of slots required for every node in Cell c to have 
transmitted successfully atleast once. Recall that p = {kiC2 log n)~^. We will bound % by stochastic 
domination. Consider a sample space S containing mutually disjoint events Ei,E2,. ■ ■ ,En^. Let 
Pr{Eg) = p{l _p)fciC2iogn-i for 1 < g < Nc. Observe that Pi > Pr{Eg) Vi G 5^ and 1 < g < Nc- 
Let E = U^i^q- We have Pe := Pr{E) = Ncp{l - p)'=iC2 iogn-i_ ^ sequence of samples be 
drawn independently from S. Let the number of samples required to be drawn from S so that each 
of the events Eg, q = 1,2, ... Nc occurs atleast once, be the random variable T^. The probability of 
occurence of -E in a given sample is Pe and hence, the waiting time in terms of number of samples 
drawn, for the event E to occur, as well as the waiting time between consecutive occurences of 
E, is given by the geometrically distributed random variable Geom(Pj;). Now, consider the events 
of successful occurences of event E. If (I — 1) distinct events among Eg, 1 < q < Nc have already 
occured, then the probability that the next occurence of E is due to an as yet unoccured event 
Eg' is (1 — ^^), as each Eg,l < q < Nc is equally probable. The number of occurences of event 
E to wait for the occurence of an as yet unoccured event among Eg,l < q < Nc is distributed 
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as Geom(l — j^)- The random variable can thus, be expressed as: = Ylf=i^'c,j where 
t;^. ~ Geom(PB) and R'^ ~ YJti Geom(l - ^). 



Now compare the following two events: (1) Event A defined as the successful transmission from 
Cell c resulting from a successful transmission by Node i in Cell c and (2) Event B defined as the 
occurrence of in a sample drawn from S due to the occurence of Eq. Observe that Pr {A) > Pr (B) . 
From this comparison, we see that 7^ will be stochastically dominated by i.e. Pr(7^ > 2) < 
Pi'(T'c — ^) G N. Further, will be stochastically dominated by the random variable Tc = 
Ylf=i ^c,j, where tcj ~ Geom(ps) and Rc ~ J2]^i Geom(l — ^) with m = [c2logn] which is an 
upper bound on Nc from (IT]). We therefore, have 



It is convenient to work with the random variable Tc because it is independent of the parameters 
of Cell c. We will obtain the moment generating functions (mgf ) of the distributions of the integer- 
valued random variables involved. Let the mgf of each random variable be denoted by the same 
character in sans serif font. For a random variable F, f{z) = Yljez P""!-^ = j) The region of 
convergence of the mgf is specified in parentheses. 



Pr {% > z) < Pr (Tc >z) Vz G N 




Psz 



S{z) {\z\>l-ps) 



l-{l-ps)z- 




5]Pr(i?, = r) [Siz)Y 





W=i {m[z-{l-ps)]-{l-l)ps) 




Thus Efp^^'^l - 



for s < log 




. Choose si = log 




. After 



some algebra, we can show the following. 
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Here c„ 



/ 2m ^ 



1 as m ^ oo by the Stirling approximation. From the Chernoff bound 
1 — ^)^^ • By the union bound, we have 



/vrm 



PS \ 
2m) 



we get Pr {% > Vi) < Pr (T^ >Vi)<c 

Pr ^maxT^ > Vi^ < MnCmyJ i^m ^1 — 
To achieve Pr(maxcgc'?c > ^i) < ^5 it is sufficient to have (1 — ^)^^ < ^^^^ 



or 



Vi > 



\ log m + log Mn + a log n — log fc + | log tt + log c„ 



-log(l 



Here, m = [c2 log n] , M„ = f^ a.rslgn l ^- Writing - log (l - = 1^ + + . . . , we can see 

that there exists a choice of V\ = O(log^n), which would be sufficient for the completion of Phase 
I, i.e., every node in every cell of the network would have successfully transmitted at least once in 
Vi slots, with probability at least (l — . 



3.2.2 Phase II: Progress to the bottom of the square 

Let the columns of cells shown in Fig.[T]be numbered Ci, C2, . . . , Q„. Let the /„ cells in each column 
be numbered from 1 to /„ from top to bottom. In this phase, we are concerned with transmissions 
in the top w := In — ^ cells of each column. In Phase I, each node has successfully received the 
transmissions by every other node in its cell. Hence, Phase II will be completed if the following 
sequence of events occurs for each column C: A successful transmission by some node in the ffist 
cell of the column, followed by a successful transmission by some node in the second cell of the 
column and so on until a successful transmission by some node in the tn-th cell of the column. 

Let the number of slots required for this sequence of events be T*-*-"^ for column C. We can see that 
T^*^) will be stochastically dominated by T'^'^'^ := YlJ=i 4*^^' "^h^re t^j^^ ~ Geom(p5'). We can thus 
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derive the following. 



T(C)(z) 



w ~w 



(l-(l-ps)^-^)"' 

{\z\ >l-ps) 
T(C), ^ Ps 

{e-^ - {1 - Ps))"" 

for s < log (^j^ 

Pr(T(^^>V2) < ^[^^ = 2-(l - ^)^^ 



rr max 

\i<j<ln J 



2 



where we have used S2 = log( \g ) in the Chernoff bound. Thus, to achieve Pr (maxi<j<;^ T^'-"^^ > V2) < 
it suffices to have (1 — < ^a!^ ow or 

a log n + log ln + w log 2 — log A; 



log(l-f ) 



V2 > 



Now, /„ = \ y 2 Tslogn "! ='"^ + 17 and hence, V2 = O slots are sufficient for the completion 

of Phase II with probability at least (l — ^) . 



3.2.3 Phase III: Progress into the sink 

Phase III comprises diffusion of the MAX into the cell containing the sink. Let the time required 
for this to happen be the random variable Tg. It is easily seen from the analysis of the sequence of 
transmission for Phase II that Pr (Tg > V3) < 2"'(1 — where w is as defined before. Calculations 
similar to those in the analysis for Phase II show that V3 = O (^y'^Iogn) ^^ots are sufficient for 
completion of this phase with probability at least (l — ^) . 



3.2.4 Bound on the overall time 

Since each of phases I, II and III get completed in O ^y'^Iogn) time slots with probability at least 

^1 — , for appropriate constants k' , the protocol One-Shot MAX achieves computation of the 

MAX at the sink in O (^y^'iogn) ^^umber of time slots with probability at least (l — ^) . If the 
protocol is followed for another V3 + V2 slots, the true MAX will diffuse to the complete bottom row. 
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and then to the complete network, the direction of diffusion being opposite to that in Phase III 
and Phase II respectively. 

3.3 Obtaining the Hop Distance 

The following algorithm Hop Distance Compute obtains the hop distance for each node in the 
network, [logd] slots are grouped into a frame and r = 0(log^ n) (r = Vi as obtained in Phase I 
analysis of protocol One-Shot MAX) frames form a superframe. The algorithm ends after (d + 1) 
superframes. 

Let the superframes be denoted by go, gi, ■ ■ ■ , gd- A node either transmits in every slot of a frame 
or it does not transmit in any slot of the frame. Each transmission is a number expressed in [log d] 
bits. At the beginning of the algorithm, the sink transmits the number expressed in [logd] 
bits in each frame of superframe go. Each node of the network other than the sink executes the 
following algorithm. Node i makes no transmission till it has decoded a transmission successfully. 
Let the first successful reception by Node i happen in a frame belonging to superframe gi and let 
the decoded transmission correspond to the number n, expressed in [logd] bits. Node i sets its 
hop distance to (nj + 1) and ignores other successfully received bits in frames from superframe g^. 
During the r frames from superframe gi+i, Node i transmits, in each frame, the number (nj + 1) 
expressed in [logd] bits, with probability p, independently of all the other transmissions in the 
network and makes no transmission with probability (1 — p). After the end of round gi+i, Node i 
makes no more transmissions. The total number of slots required is {d + l)T[log d] . 

Lemma 1 The nodes of the network correctly compute their minimum hop distance from the sink, 
using Hop Distance Compute in 0{y/n log^ n) time slots with probability at least (l — ^) for any 
positive a and some constant k. 

We omit the proof of this lemma. 

3.4 Proof of Theorem [2] 

Let the set of nodes at hop distance h be Gh- Let ti^r be the first slot in round r that Node i 
transmits succesfully in. The number of slots in a round is r = 0(log^n) (r = Vi from Phase I). 
Every node in the network would have transmitted successfully at least once in each round of r 
slots w.h.p. Let /imax < dhe the largest hop distance of a node in the network. In the proof, we 
will assume that each node of the network transmits successfully in each round at least once. We 
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claim that 



maxrj(r, ti^r) = max Zj{r — d + h) 



for < < /imax and r > d — h. The sink being at hop distance 0, proving the claim will complete 
the proof. Assume that the claim is true for ho < h < /imax for r > d — h. We shall show that the 
claim will then be true for h = ho and for r > d — ho- Consider transmissions by the nodes at hop 
distance ho in round (r + 1). 



max Ti{r + l,ti^r+i) = max {max{Zi(r + 1 — d + ho),Yi{r)}} 



Since each node at hop distance (ho + 1) transmits successfully at least once in round r, the 
transmission of each such node is decoded successfully by some node at hop distance ho- Hence, 

max l^(r) = max Tj(r,tj^r) 

= max — d + ho + 1) 

ieUh(,+i</<d'j'/ 

where the second equality follows from the induction hypothesis. Hence, 

max Tj(r + 1, = max{ max Zi{r + I ~ d + /iq), 

max Zj (r — d + /lo + 1)} 

jeUho+i</<d Gj 

= max Zj(r — d + /iQ + 1)} 

jeUft(,</<d 

which proves the claim for hop distance ho for round (r+1). By induction, the claim is true for each h 
and each round r > d—h. Therefore, the sink Node s correctly sets Z{r—d) = max{ (r— d), ys(r)}. 
The delay of the protocol is dr = 0(\/n log^ n) slots. 

As transmissions by different nodes are independent, the analysis in the diffusion of phase I of 
One-Shot MAX carries over. The probability that the computed value of Z{r) is incorrect for any 
given round is upper bounded by for any constants a, A; > 0. 
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4 Discussion 



The total number of transmissions (successful as well as unsuccessful) in one execution of One-Shot 

3 /2 

MAX is Q( j "3/2 )• In Pipelined MAX, a total of 0(nlogn) transmissions are made per round. Note 
that the corresponding number is Q(n) with a coordinated protocol for both cases. 

Our analysis can be extended to the case where the nodes use pure Aloha as the MAC. We need 
to use a transmission rate rather than a transmission probability. The success probabilities are 
calculated similarly except that we now have a collision window that is twice the packet length. 
All calculations are analogous. 

It is fairly straightforward to show that in a noiseless, structure-free broadcast network, the his- 
togram can be computed in 0(n) slots w.h.p. In the noisy broadcast network, by a simple modi- 
fication of the protocol of [Tj, we can show that the histogram can be computed in O(nloglogn) 
slots w.h.p. 
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